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ABSTRACT 

It was recently proposed that metal-rich white dwarfs (WDs) accrete their metals from compact debris 
disks found to exist around more than a dozen of them. At the same time, elemental abundances 
measured in atmospheres of some WDs imply vigorous metal accretion at rates up to 10 11 g s _1 , 
far in excess of what can be supplied solely by Poynting-Robertson drag acting on such debris disks. 
To explain this observation we propose a model, in which rapid transport of metals from the disk 
onto the WD naturally results from interaction between this particulate disk and spatially coexisting 
disk of metallic gas. The latter is fed by evaporation of debris particles at the sublimation radius 
located at several tens of WD radii. Because of pressure support gaseous disk orbits WD slower than 
particulate disk. Resultant azimuthal drift between them at speed < 1 m s _1 causes aerodynamic 
drag on the disk of solids and drives inward migration of its constituent particles. Upon reaching the 
sublimation radius particles evaporate, enhancing the density of metallic gaseous disk and leading to 
positive feedback. Under favorable circumstances (low viscosity in the disk of metallic gas and efficient 
aerodynamic coupling between the disks) system evolves in a runaway fashion, destroying debris disk 

on time scale of ~ 10 5 yr, and giving rise to high metal accretion rates up to Mz ~ 10 10 - 
in agreement with observations. 

Subject headings: White dwarfs — Accretion, accretion disks — Protoplanetary disks 
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1. INTRODUCTION. 

Recent detections of near-infrared excesses around a 
number of metal-rich white dwarfs (WD) imply existence 
of warm circumstellar material reprocessing stellar radi- 
ation (Zuckerman & Becklin 1987; Graham et al. 1990; 
Farihi et al. 2010). Spectral modeling suggests that this 
material resides in an extended, compact, optically thick 
and geometrically thin disk (Jura 2003; Jura et al. 2007), 
similar to the Saturn's rings (Cuzzi et al. 2010). Disks 
have inner edges at several tens of the WD radii i?*, 
roughly consistent with being set by particle sublima- 
tion at these locations. Their outer radii lie close to 
the Roche radius of the WD Rr ~ 1 R Q , supporting the 
suggestion by Jura (2003) that such compact debris disks 
are produced by tidal disruption of asteroid-like bodies 
scattered into low-periastron orbits by gravitational per- 
turbations of massive planets, which have survived the 
AGB phase of stellar evolution. 

Availability of large reservoir of high-Z elements in the 
form of debris disk in the immediate vicinity of some 
WDs naturally led to the suggestion (Jura 2003) that 
metal enrichment of these stars is caused by accretion 
from such disks. This scenario provides a promising al- 
ternative to the previously proposed interstellar accre- 
tion model of metal pollution of WDs (Dupuis et al. 
1993), which is not consistent with observations of WDs 
with He atmospheres. 

Theoretical estimates of settling time of heavy ele- 
ments in WDs imply that their observed atmospheric 
abundances can be maintained against gravitational set- 
tling by accretion of metals at rates Mz ~ 10 6 — 10 11 
g s _1 (Farihi et al. 2009, 2010). If the circumstellar ac- 
cretion hypothesis is correct, an evolving disk of debris 
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must be able to supply such Mz to the WD. 

The actual transfer of metals from the disk of solids 
truncated at sublimation radius R s to the WD must 
be accomplished in this picture via the gas disk extend- 
ing from the WD surface to R s and beyond. Observa- 
tional evidence of such gaseous component around sev- 
eral metal-rich WDs hosting compact debris disks has 
been found by Gansicke et al. (2006, 2007, 2008) in the 
form of double-peaked emission lines of Ca II and Fe II. 
These spectroscopic signatures are naturally explained 
as arising in a disk of metallic gas (no H or He emission 
lines have been detected around these WDs) in Keple- 
rian rotation around WDs and spatially coincident with 
dusty disks (Melis et al. 2010). 

Even though metals are passed to the WD through 
the gaseous disk the rate of mass transfer Mz is ulti- 
mately controlled by evolution of the particulate disk. 
By analogy with Saturn's rings one expects that the col- 
lisional viscosity in debris disk is too low to drive the 
non- negligible Mz- However, radiation of the WD can 
be quite important and Rafikov (2011; hereafter Rll) 
demonstrated that Poynting-Robertson (PR) drag on the 
disk naturally drives mass accretion at the rate Mpr > 
10 8 g s" 1 . While this M z is high it still falls short of 
explaining the highest observed Mz ~ 10 10 — 10 11 g s _1 . 

Here we propose a new picture of the debris disk evo- 
lution, which naturally combines several physical ingre- 
dients present in close vicinity of disk-bearing WDs. 

2. DESCRIPTION OF THE MODEL. 

Our model of the WD-disk system is described below 
and is schematically illustrated in Figure [1] We assume 
a disk of solid particles to extend from the Roche radius 
Rr all the way in to the sublimation radius 
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where T s is the sublimation temperature of solid particles 
(which we take to be ~ 1500 K if particles are silicate), 
T*.4 = T*/(10 4 K) is the normalized stellar temperature 
T*, R± is the WD radius, and e is the ratio of particle 
emissivities for starlight and for its own thermal radiation 
(in the following we assume macroscopic particles and set 
e = 1). For R+ sa O.Oli?© (Ehrenreich et al. 2011) one 
finds R s ~ 0.2 Rq, in agreement with observations (Jura 
et al. 2007, 2009). 

At present the size of particles a constituting the disk 
is rather poorly constrained. Graham et al. (1990) sug- 
gested a s» 10— 100 cm. At the same time high-resolution 
infrared spectroscopy with IRS onboard of Spitzer reveals 
strong 10 /im emission feature in disk spectra indicative 
of the existence of a population of small, micron-sized sil- 
icate dust particles (Jura et al. 2009). But the fraction of 
total disk mass locked up in such fine dust is unknown. 
The actual particle size is not very important for our 
treatment as long as the disk is optically thick. Then it 
is similar to Saturn's rings and most likely behaves as a 
granular flow. In the following we will treat particulate 
disk as if it were a solid plate. 

Solid particles brought to R s sublimate feeding a disk 
of metallic gas at this location. Viscous torques cause 
it to spread all the way to the WD surface providing 
means of metal transport from R s to the star. However, 
because of angular momentum conservation part of the 
newly produced metallic gas has to move outwards of R s 
(Lynden-Bell & Pringle 1974) naturally explaining the 
existence of gaseous disks spatially coexisting with dusty 
debris disks in some systems (Melis et al. 2010). Jura 
(2008) proposed collisional sputtering of small asteroids 
as another source of metallic gas; for the sake of clarity 
we do not consider this mechanism here. 

We emphasize here that the two disks have an over- 
lap in radial distance, while in vertical direction gaseous 
disk is much more extended than the disk of particles, see 
Figure [TJ The latter must be very thin because inelas- 
tic collisions between particles rapidly damp any vertical 
random motions. 

External gas disk outside R s has temperature different 
from (higher than) that of the dust disk at the same 
radii (Melis et al. 2010) because of the different balance 
of heating and cooling for the two disks. This keeps gas 
temperature T g above the sublimation point even outside 
R s , although some condensation of metallic gas on the 
particle surfaces may be happening there (we neglect it 
in this work). 

Simultaneous existence of the two disks of high-Z ele- 
ments in different phases drives their mutual evolution in 
the following way. Any coupling between the outer por- 
tion of the gas disk and particulate disk acts to transfer 
angular momentum from the faster rotating particulate 
disk to the slower spinning gaseous disk, causing inward 
motion of particles in the disk of solids. If the coupling 
is strong enough positive feedback becomes possible in 
the system: increasing mass of gaseous disk leads to the 
increase of Mz through the disk of solids (as described 
below), which in turn reinforces evaporation at R s and 
increases gaseous mass even further. 

Coupling between the gaseous and particulate high-Z 
disks arises because gaseous disk orbits WD at angular 
speed fig slightly lower than the Keplerian speed £Ik at 




Fig. 1. — Schematic picture of two interacting disks (black - 
particulate, gray - gaseous) around the WD. Disk of solid parti- 
cles is confined between the sublimation radius R s (where solids 
evaporate feeding gas disk) and Roche radius Rr. Flow of high-Z 
material is indicated by arrows (white - solid particles, black - gas). 
Inward migration of solids is induced by drag exerted on them by 
the gas disk, which orbits WD slower than the disk of particles. 



which the disk of solid particles rotates. This difference 
is caused by the pressure support present in gaseous disk, 
Vtg—Q,K ~ (2flKrp)~ 1 dP/dr, and is known to cause a va- 
riety of important effects in protoplanetary disks, such as 
inward migration of solids (Weidenschilling 1977). Rela- 
tive azimuthal velocity between the gaseous and partic- 
ulate disks at distance r from the WD is 

where rj ~ 1 is a constant, M*^ = M+/Mq is the nor- 
malized WD mass M*, /j,2s is the mean molecular weight 
of the metallic gas normalized by 28m p (value of /i for 
pure Si), T g ^ = T g /(10 3 K) is the normalized gas tem- 
perature T g: and c s ~ 0.5 km s _1 (Tg i3 //Z28) 1 ^ 2 is the gas 
sound speed (clearly v v -C c s ). 

This azimuthal drift gives rise to aerodynamic drag 
between the disks. Azimuthal drag force f v per unit sur- 
face area of a particulate disk with surface mass density 
Tid causes inward radial migration of the disk material at 
speed v r = 2f v / (f2jfE<j). This inward particle drift gives 
rise to mass transport at the rate 

M Z = 27rnvE d - (3) 

Note that if f v is independent of then the same is 
true for Mz- 

External force f v per unit area can be generally repre- 
sented in the form 

f v = AUg + B, (4) 

where the first term describes the coupling between the 
gas disk with surface density E g and particulate disk, 
and constant A determines the strength of coupling. It 
is natural to expect that drag scales with gas density E s ; 
this expectation is confirmed in 2J 

Second term B represents forces acting even in the ab- 
sence of coupling to gas disk. PR drag is an example of 
such force and one can easily show (RH) that for PR 
drag Bpr s» a^KL^(f> r /(4:Trrc 2 ), where c is the speed of 
light, <j) r ~ 1 is the efficiency of radiative momentum 
absorption by the disk surface, and a ps (4/Sw)(R is /r) 
(Friedjung 1985) is the incidence angle of stellar radia- 
tion. From equation ([3]) the rate of mass transport due 
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to PR drag alone (when A = 0) at R s is (Rll) 



M PR = - 



ao 8 g s~\ 
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(5) 
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In the following we will assume that B = Bp a. 

3. COUPLED EVOLUTION OF THE PARTICULATE AND 
GASEOUS DISKS. 

Particle sublimation at R s increases the mass of 
gaseous disk at the rate Mz- Assuming that surface 
density of gaseous disk S g varies on scale ~ R s we 
can write that sublimation increases £ g at the rate 

£+ ~ M/(ttR 2 s ) = 4f v /(n K R s ). At the same time, 

viscous spreading reduces £ g at the rate £_ ~ Ti g /t u , 
where t v is the characteristic viscous time in the disk 
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assuming a-parametrization of viscosity v = aCg/ilx 
(Shakura & Sunyaev 1973). This timescale can be very 
short if a is not very small. 

We can now describe the evolution of £ g in the vicinity 
of R s with the following heuristic equation: 



dt 
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A solution of this equation satisfying the initial condition 
T, g (t = 0) = (no gas disk initially) is 



s fl (t) = 
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where we used equation ([5]) and 
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is the sublimation time on which the mass of 

gaseous disk increases if Bpp = 0. From equations ([3]), 
((H) and © the rate at which particles sublimate is 
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and the total mass lost by the debris disk to sublimation 
M s can be easily obtained by integrating this expression. 

The behavior of Mz and M s is shown in Figure [5] for 
different values of t s /t v , which clearly demonstrates that 
for t < t s (i.e. as long as M s < Mppt s ) debris disk 
evolution is insensitive to t s /t v and Mz ~ Mpp. This 
is because initially the disk of metallic gas is not dense 
enough for the drag it produces on the debris disk to 
compete with the PR drag — it takes certain time to 
accumulate enough gas mass by particle sublimation. 

This means, in particular, that if the debris disk starts 
with mass M^, which is lower than the critical mass 
Mppt s , then its evolution is determined only by the PR 
drag and coupling to the gaseous disk is never effective — 




Fig. 2. — Time evolution of (a) Mz and (b) total sublimated 
mass M s . Evolutionary tracks are shown for different values of the 
feedback parameter t^/t s = T and demonstrate runaway behavior 
for tujts > 1 and saturation of M% for t^/ts < 1. Note that for 
t < t s the behavior of both Mz and M s is independent of tt//t s . 



the debris disk is eroded before gas disks grows massive 
enough. Lifetime of the debris disk is then 
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independent of the relation between t s and t v . 

In the opposite case of a massive initial disk — Md 
Mppt s — evolution does depend on t s /t v . Whenever 
tv J; the action of viscosity is so effective at remov- 
ing gas released by particle sublimation at R s that the 
gas does not accumulate there. Then E g simply satu- 
rates at the constant low level, accretion rate tends to 
M z w M PR (1 - t v jt s y x - M PR for t > t v , and disk 
gets exhausted on timescale tpp given by equation (fTTj) . 
This evolutionary path clearly cannot explain the highest 
observed values of Mz- 

However, in the case of a massive disk and t v > t s 
sublimation supplies gas to the r ~ R s region faster than 
viscous diffusion removes it, and E ff grows exponentially 
on timescale ~ t s (as long as particle disk has enough 
mass to provide the source). This means that Mz also 
increases exponentially for t > t s (see Figure [2|) . This 
runaway stops and gaseous disk wanes only when the 
mass of particulate disk Md is exhausted. The latter 
happens on the runaway timescale 



tr 
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determined by the condition M s (t mn ) 
M d . Clearly, t run < t PR if M d /{M PR t s 



= J*' un M z dt 
> 1. 
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Thus, the ratio 
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t. 



(13) 



is the critical parameter determining the strength of pos- 
itive feedback for massive disks: if J- > 1 a particulate 
disk gets rather rapidly (within several t s ) converted into 
metallic gas at R s , with its subsequent viscous accretion 
onto the WD on timescale of several t v . Whereas for 
T < 1 feedback is not strong enough to reinforce gas 
supply at R s and disk slowly evolves on timescale tp R 
due to PR drag. As equation (fT3")) shows, runaway evo- 
lution and high Mz require (1) strong coupling between 
the gaseous and particulate disks and (2) low viscosity. 

4. AERODYNAMIC COUPLING. 

Relative azimuthal motion between gaseous and par- 
ticulate disks at speed v v induces aerodynamic drag be- 
tween them. Drag produces force per unit area of the 
disk (Schlichting 1979) 



(14) 



where p g = Q.K^ g /c s is the gas density and Re" 1 is a 
proportionality constant. There is a significant spread of 
opinions regarding the value of Re* characterizing drag 
by the turbulent flow on a smooth solid plate, with num- 
bers ranging between Re* ss 20 (Dobrovolskis et al. 1999) 
to Re* w 500 (Goldreich & Ward 1973). It is also likely 
that a smooth plate approximation underestimates drag 
(overestimates the value of Re*) for the particulate debris 
disk, which does not have continuous surface and may in- 
teract with gas more like a rough plate (Schlichting 1979) 
or even as a combination of individual particles, in which 
case smaller Re* is more appropriate. 

According to equation (|14j) aerodynamic drag force can 
be written in the form / = AT, g with 
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where equation @ was used. From equation ^ subli- 
mation time at r — R s is 



Re* GM* 
4ry 2 c 3 
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where c s ,i = c s /(l km s 1 ). The critical mass separating 
low- and high-mass disks (see S|3|) is 
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The feedback parameter for aerodynamic coupling is 
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and depends on both the viscosity parameter a and the 
strength of aerodynamic coupling, parameterized by Re*. 

Equation ([T"8|) shows that runaway evolution with T > 
1 requires rather low viscosity in the gaseous disk, at the 
level of a ~ 10~ 3 — 10 -4 , and for higher Re* smaller 
a is needed. Viscosity is most likely provided by the 
magnetorotational instability (MRI) , which requires cer- 
tain level of ionization to operate effectively. In the ideal 



MHD limit simulations with no net flux typically produce 
a ~ 10~ 2 (Hawley et al. 1995). However, in our scenario 
gaseous disk exists in immediate contact with the partic- 
ulate disk, which is observationally known to contain a 
population of micron size dust grains (Jura et al. 2009). 
MRI-driven turbulence will mix some of this fine dust 
with the gas, lowering ionization fraction (small grains 
have large surface area and are very efficient charge ab- 
sorbers), and giving rise to the non-ideal MHD effects, 
e.g via the increased resistivity (Balbus 2009). The lat- 
ter are known from simulations to decrease effective a 
substantially, down to a ~ 10~ 4 or lower (Fleming et 
al. 2000; Bai & Stone, in preparation). As a result, it is 
conceivable that a in our model can be much lower than 
in the ideal MHD limit of MRI. 

Thus, presence of dusty debris disk in close contact 
with the gaseous disk can quite naturally lengthen vis- 
cous timescale t v and facilitate gas accumulation, mak- 
ing possible runaway evolution due to aerodynamic drag. 
We note that after the particulate disk completely sub- 
limates the value of a in the resultant gaseous disk will 
go up since the source of small dust particles lowering 
ionization has disappeared. Closer to the WD, in the 
dust- free region r < R s , ionization and a can always be 
higher than in the external disk. Faster viscous evolu- 
tion in this region (resulting in lower E g ) may explain 
lack of line emission from the inner part of the gas disk 
(Gansicke et al. 2006, 2007, 2008; Melis et al. 2010). 

5. DISCUSSION. 

Our calculations clearly demonstrate that compact de- 
bris disks around WDs are self-destructive whenever 
there is strong coupling between the gaseous and solid 
components. Adopting for illustration Re* = 20 and set- 
ting all other dimensionless constants to unity we find 
from equations (|l6]) - (pf t s « 2 x 10 4 yr, M PR t s « 
6 x 10 19 g and that a < 10~ 4 is needed for T > 1. 

Then the typical lifetime of a massive disk with Md = 
10 22 g (mass of a 200 km asteroid) in the runaway sce- 
nario (if the viscosity is low) is on the order of several 
sublimation timescales, i.e. about 10 5 yr. According to 
equations ()10j) and (|12p in this case the maximum Mz for 
this Md, achieved right before the debris disk completely 
disappears, is max(Mz) ~ Md/t s ~ 10 10 g s _1 . These 
numbers agree with observational inferences (Farihi et al. 
2009, 2010) quite well. 

However, if positive feedback is not strong enough to 
drive runaway (e.g. because a > 10~ 4 or higher Re*), 
then Mz ~ Mpr and it takes 3 Myr for the same disk 
to be exhausted, as equation (fTTj) demonstrates. In this 
case gas surface density at R s saturates at the low level 
E g ~ M PR t s / [ttR 2 s ) ~ 0.1 g cm -2 , see equation ©. 

Equation (fTTj) also shows that a low-mass disk with 
Md = 10 19 g (mass of a 20 km asteroid) is destroyed 
by PR drag very rapidly, within several thousand years 
(as long as the disk is optically thick to incoming stellar 
radiation, see Rll). 

Our model of debris disk evolution naturally explains 
coexisting gaseous and particulate debris disks reported 
in Melis et al. (2010). Systems with reported IR excesses 
but lacking emission lines of high-Z elements in gas phase 
may simply possess lower mass gaseous disks, which have 
not had enough time to develop by sublimation of solids 
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at R s . Compositional variations between different WDs 
may also explain such systems. 

Calculations presented in this work are rather sim- 
ple and local in nature. We studied only one coupling 
mechanism — aerodynamic drag, while other possibil- 
ities may be available as well (e.g. induction interac- 
tion (Drell et al. 1965; Gurevich et al. 1978) between the 
MRI-generated B-field in the gaseous disk and the de- 
bris disk particles). Also, here we did not consider low 
surface density disks (we use only the thin plate approx- 
imation), non-trivial initial radial distribution of debris 
surface density, fate of angular momentum lost by debris 
disk and absorbed by gas disk, and so on. Future global 



models of coupled evolution of gaseous and particulate 
debris disks (Bochkarev & Rafikov, in preparation) will 
take these issues into account to provide a more accurate 
description of WD pollution with circumstellar high-Z 
material. 
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